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ABSTRACT 

Aims. Building on the two-point correlation function analyses of the VIMOS Public Extragalactic Redshift Survey (VIPERS), we 
investigate the higher-order correlation properties of the same galaxy samples to test the hierarchical scaling hypothesis at z ~ 1 and 
the dependence on galaxy luminosity, stellar mass, and redshift. With this work we also aim to assess possible deviations from the 
linearity of galaxy bias independently from a previously performed analysis of our survey. 

Methods. We have measured the count probability distribution function in spherical cells of varying radii (3 < R < 10/7 1 Mpc), 
deriving crg 9 (the galaxy rms at 8 h~ l Mpc), the volume-averaged two-, three-, and four-point correlation functions and the normalized 
skewness S 2g and kurtosis S^ g for different volume-limited subsamples, covering the following ranges: -19.5 < M B (z = 1.1) - 
5 log(/r) < -21.0 in absolute magnitude, 9.0 < log (M,/M e h~ 2 ) < 11.0 in stellar mass, and 0.5 < z < 1.1 in redshift. 

Results. We have performed the first measurement of high-order correlation functions at z ~ 1 in a spectroscopic redshift survey. Our 
main results are the following. 1) The hierarchical scaling between the volume-averaged two- and three-point and two- and four-point 
correlation functions holds throughout the whole range of scale and redshift we could test. 2) We do not find a significant dependence 
of S 3g on luminosity (below z = 0.9 the value of S 2g decreases with luminosity, but only at lcr-level). 3) We do not detect a significant 
dependence of S 2g and Sn g on scale, except beyond z ~ 0.9, where S 2g and Sn g have higher values on large scales (R > IO/7 1 Mpc): 
this increase is mainly due to one of the two CFHTLS Wide Fields observed by VIPERS and can be explained as a consequence 
of sample variance, consistently with our analysis of mock catalogs. 4) We do not detect a significant evolution of S 2g and S 4 g with 
redshift (apart from the increase of their values with scale in the last redshift bin). 5) a% g increases with luminosity, but does not 
show significant evolution with redshift. As a consequence, the linear bias factor b = crg g /crg m , where erg,,, is the rms of matter at a 
scale of 8 hr 1 Mpc, increases with redshift, in agreement with the independent analysis of VIPERS and of other surveys such as the 
VIMOS-VLT Deep Survey (VVDS). We measure the lowest bias b = 1.47 ± 0.18 for galaxies with M B (z = 1.1) - 51og (h) < -19.5 
in the first redshift bin (0.5 < z < 0.7) and the highest bias b = 2.12 ± 0.28 for galaxies with M B (z = 1.1) - 51og(/i) < -21.0 in the 
last redshift bin (0.9 < z < 1.1). 6) We quantify deviations from the linear bias by means of the Taylor expansion parameter b 2 . We 
obtain b 2 = -0.20 ± 0.49 for 0.5 < z < 0.7 and b 2 = -0.24 ± 0.35 for 0.7 < z < 0.9, while for the redshift range 0.9 < z < 1.1 we find 
b 2 = +0.78 ± 0.82. These results are compatible with a null non-linear bias term, but taking into account another analysis for VIPERS 
and the analysis of other surveys, we argue that there is evidence for a small but non-zero non-linear bias term. 
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1. Introduction 

In the standard model of structure formation, the growth of 
density fluctuations from a primordial Gaussian density field is 
driven by gravity; it is possible to follow the evolution of these 
fluctuations through analytical and numerical approaches and 
predict the statistical properties for the dark matter field and dark 
matter haloes. Galaxies form in a complex process following the 
baryonic infall into dark matter halos: this means that the com¬ 
parison between theory and observations is not straightforward, 
but it also implies that the spatial distribution of galaxies con¬ 
tains a wealth of information relevant for both cosmology and 
the physics of galaxy formation. 
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>From the analysis of the VIMOS-VL T Deep Survey, based 
on spectroscopic rcdshift s.lMarinoni et al. (i2005l) (see also Mari- 
noni et al. 120081 ) found that the value of S 3 for luminous (M# < 
- 21 ) galaxies is consistent with the local value at z < 1 while 
decreasing beyond z ~ 1 , and that the bias is non-linear. 

In this paper we analyse the high-order correlations and mo¬ 
ments of the first release of the VIMOS Public Extragalactic 
Redshift Survey (VIPERS0) in the redshift range 0.5 < z < 1.1 
as a function of luminosity and stellar mass. We also derive an 
estimate of the non-linear bias. Our analysis extends those pre¬ 
sented in a number of recent works that have investigated various 
aspects of galaxy clustering in the VIPERS sample. Some works 
have focused on two-point statistics, like the standard galaxy- 
galaxy two -point correlation functio n to estimate redshift space 
distortions dde la Torre et al.i l2013h and its evolution and de¬ 
pendence on galaxy properties (iMarulli et al.ll20l3h . A different 
type of t wo-point statistics, the clustering ratio, has been intro¬ 
duc ed bvl Bel & M arinonil (1201 4h and applied to VIPERS galax¬ 
ies ( Bel et aljt'20 1 -417 to estimate the mass density parameter Q w . 
iMicheletti et al.i ( 20141) have searched the VIPERS survey for 
galaxy voids and characterized their properties by means of the 
galaxy-void cross-correlation. Bel et al. (2015, in preparation) 
have proposed a method to infer the one-point galaxy^ probabil¬ 
ity function from counts in cells that lDi Porto et al .i (12014(1 have 
exploited to search for and detect deviations from linear bias; a 
resul t that we directly compare our results with. Finally, Cucciati 
et al. (i2014li *studied different methods for accounting for gaps in 
the VIPERS survey and assessing their impact on galaxy counts. 

As cosmological parameters we have adopted Hq = 70 
km/s/Mpc, CIm — 0.25, Qa = 0.75, but all cosmology- 
dependent quantities are given in Hq = 100 km/s/Mpc units as¬ 
sociated with the corresponding power of h = /-/ ( , / 100 . 


2. High-order statistics 

In this section we resume the formalism and define the statistical 
quantities measured in our work. 

The volume-averaged /-point correlation functions are 
given by 

h( V) = ^jffj dv i- dv J> ( 2 ) 

where for spherical cells (used in this work) £, is a function 
of the cell radius R and V = AttR '/?>. 

The volume-averaged two-point correlation function gives 
the variance of the density contrast: 


The volume-averaged correlation functions correspond to 
the reduced moments and up to the fourth order are given by 
the following relations: 

1 

£2 = H2 - = 

N 

£3 -F3 -3= + — 

N /W 

£t =m -6— + 11^7 - 3jU, - —r. (5) 

N N N 

An alternative way to estimate the high-order correlations is 
through the factorial moments m k \ 

m k = P(N)NK (6) 

N 

where 


N k = N(N- l)...(N-k + 1) (7) 

_is the falling factorial power of order k (see e.g. lGrahamet al.i 

11994 . 

In fact, for a local Poisson process the moments about the 
origin of a stochastic field are given by the factorial moments of 
/V; as our variable is the number density contrast {N — N)/N, we 
have to convert the factorial moments m k into the moments about 
the mean (central moments) /j' k through the standard relations 

/4 = m 2 ~ m\ 

yL /3 = m 3 — 3wiW2 + 2 m j 

yL/4 = »4 — 4;HpH3 + 6m\m.2 - 2>m\. (8) 

We can finally derive the volume-averaged /-point correla¬ 
tion functions 


T-^J 
w - — 

N 


and the normalized moments S , 



( 9 ) 


( 10 ) 


The normali zed moments can also be obtained through a re - 
cursive formula dSzanudi & Szalav)ll993 l lColombi et aklboOdl) : 


a 2 (R) = &(*)■ 


(3) 


The volume-averaged /-point correlation functions can be 
easily derived from the moments of the count PDF P(N, R), that 
is, the probability of counting N obj ects in a randomly chosen 
spherical volume of radius R (see (Peebles! 198(j) . For simplicity, 
in the following we omit the dependence on R. At a fixed scale 
R, the centred moments of order / are 


RJ 


J] P(N ) 

N 


' n-n Y 

. N , 


(4) 


where N = nV = £ NP(N) is the mean number of objects in a 
cell of radius R. 

2 http://vipers.inaf.it 


% 2 N j - 1^ /! (/ - k)S j-k'n k 

n j c J h (J ~ k) - k ■ N c 

where 


( 11 ) 


N c = Nf 2 . (12) 

The values given in this paper were calculated using factorial 
moments. 

At a fixed scale R, the deterministic bias parameter b can 
be directly measured through the square root of the ratio of the 
galaxy variance cr 2 to the matter variance cr 2 n : 


b(z) = 


°~g(z) 
O" m(z) 


(13) 
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In the case of linear biasing, the galaxy density contrast 5 g 
is proportional to the matter density contrast 6 m by a constant 
factor b, 5 g = b6 m : there is no dependence on scale, and b is the 
only parameter that completely defines the relation between the 
galaxy and matter distribution. 

As we have noted in the introduction, the linear biasing can¬ 
not have a general validity. It is more general and realistic to 
assume a local, deterministic non-li near bias b(z, 6 m , R), wh ich 
can be written as a Taylor expansion dFrv & Gaztanagalll993h : 


e V bk Jc 

5 < = Lp 5 - 


k =0 


(14) 


where b\ = Zt. lFrv & Gaztanagalri 19931) have shown that such 
a local bias transformation preserves the hierarchical properties 
of the underlying matter distribution in the limit of small fluctu¬ 
ations (large scales). 

In the case of linear bias, bk = 0 for A: > 1, and the galaxy and 
matter normalized moments are then related by the following 
equation: 


s„ = f&. 05) 

In general, the deviation from linear biasing is measured by 
taking the second order of the expansion. In this case, the galaxy 
normalized skewness is given by the following relation: 


S 





(16) 


3. VIPERS survey 

The VIMOS Public Extragalactic Redshift Survey (VIPERS) is 
an ongoing ESO Large Programme aimed at determining red- 
shifts for ~ 10 5 galaxies in the redshift range 0.5 < z < 1.2, to 
accurately and robustly measure clustering, the growth of struc¬ 
ture (through redshift-space distortions) and galaxy properties 
at an epoch when the Universe w as about half its c urrent age 
JGuzzo & The Vipers Teamll2013 l: lGuzzo et al.ll2014li . The sur¬ 
vey is divided into two separate areas and will cover ~ 24 deg 2 
when completed. The two areas are the so-called W 1 and W4 
fields of the Canada-France-Hawaii Telescope Legacy Survey 
Wide (CFHTLS-Wide); the CFHTLS optical photometric cat¬ 
alogue^ constitute the parent catalogue from which VIPERS 
spectroscopic targets were selected. The VIPERS survey strat¬ 
egy is optimiz ed to achieve a good c ompleteness in the largest 
possible area (tScodeggio et al.1120091) . Galaxies are selected to 
a limit of Gb < 22.5, further applying a simple and robust gri 
colour pre-selection to effectively remove galaxies at z < 0.5. 
In this way, only one pass per field is required, allowing us to 
double the galaxy sampling rate in the redshift range of interest 
with respect to a pure magnitude-limited sample (~ 40%). The 
final volume of the survey will be 5 x 10 7 /z 3 Mpc 3 , comparable 
to that of the 2dFGRS at z ~ 0.1. 

VIPERS spectra are obtained using the VLT Visib le Multi- 
Object Spectrograph (VIMOS. iLe Fevre et al.ll2002l Le Fevre 
et al. l2003h at moderate resolution (R — 2101. with the LR Red 
grism at R = 210 and a wavelength coverage of 5500-9500A. 

3 Mellier, Y., Bertin, E., Hudelot, P., et al. 2008, 
http://terapix.iap.fr/cplt/oldSite/Descart/CFHTLS-T0005-Release.pdf. 


The typical radial velocity error is 140(1 + z) km sec -1 . A dis¬ 
cussion of the survey data redu ctio n and the first management 
infrastructure were presented in iGarilli et al.l (12012l) a nd the de¬ 
tailed description of the survey was given bv lGuzzo et al.l (l2014lf . 

The data set used in this and the other published papers is 
the VIPERS Public Data Relea se 1 (P DR-1 ) cata logue, made 
available to the public in 2013 (IGarilli et al.l 20141) . It includes 
about 47,000 reliable spectroscopic redshifts of galaxies and ac¬ 
tive galactic nuclei (AGNs). We here only selected galaxies with 
reliable redshift, that is, with spectroscopic quality flags 2,3,4, 
or 9 (see IGarilli et al.ll2014l for the definition). 

To avoid regions dominated by large gaps, we here selected 
a subset of the total area covered by VIPERS: our limits are 
02 /! 01"'00' 5 < RA < 02*34 m 50 s , -5.08° < DEC < -4.17° 
(7.67 square degrees) in W 1 and 22 h 0V”l2 s < RA < 22 /, 18 m 00 I , 
0.865° < DEC < 2.20° (5.60 square degrees) in W4. 

We defined volume-limited subsamples with different abso¬ 
lute magnitude and stella r mas s limits, following the same cri¬ 
teria as in iMarulli et al. 1 ( 20131) . The choice of these particular 
samples is discussed in detail in that paper; here we recall their 
main properties. 

The rest-frame B-band absolute magnitude and the stellar 
mass were estimated through the HYPERZMASS program (Bol- 
zonella et al. 200(1 Bolzonella et al.l2010t) . which applies a spec¬ 
tral energy distribution (SED) fitting technique. To take into ac¬ 
count luminosity evolution, we fixed as a reference limit the lu¬ 
minosity at our maximum redshift (z — 1.1) and ass umed an 
evolution M(z ) = M( 0) - z (see Meneux et all 120091 and also 
lllbert et alJl2005llz ucca et 

We did not correct the mass limit of the stellar-mass-limited 
subsamples; this limit was therefore kept fixed within each red¬ 
shift b in because the evolution of M, is negligib l e in our redshift 
range dPozzetti et al.1120071 iPozzetti et al.lUOlOl iDavidzon et al.l 
12013 ) 


The respective numbers of galaxies for the different subsam¬ 
ples are given in TableQ] We note that these numbers are slightly 
different from those in lMarulli et al.l d2013l) because we applied 
more stringent angular limits to avoid regions nearby promi¬ 
nent gaps that might affect the counts in spherical cells (while 
the direct estimate of the two-point correlation function through 
counts of galaxy pairs can be easily corrected for by using a ran¬ 
dom catalogue with the same survey geometry). 


4. Analysis of mock catalogues 

We used mock catalogues derived from cosmological simula¬ 
tions to estimate not only the statistical errors and the uncer¬ 
tainty related to cosmic variance, but also the systematic errors 
that are due to the inhomogeneous spectroscopic completeness 
and the specific geometry of the two fields. A detailed descrip¬ 
tion of the wav these mocks were built was given by dc la Torre 
et al. (120131) . 

We analysed a set of 26 independent mock catalogues based 
on the dar k matt er halo catalogue of the MultiDark simulation 
(iPrada et al.l2012l) . which assumes a flat AC DM cosmology with 
(Dm, O a , Q. b , h, n, cr 8m ) = (0.27, 0.73, 0.0469, 0.7, 0.95, 0.82). 
This catalogue was populated with galaxies using hal o occupa¬ 
tion di stribution prescriptions, as described in Ide la Torre et al.l 
(l2013i) . In particular, the original halo catalogue was repopulated 
with halos_below the resolution limit with the new technique of 
Ide la Torre & Peacockl (12013 ). which enables reproducing the 
range in stellar mass and luminosity probed by VIPERS data. For 
luminosity-limited subsamples, galaxy luminosities were cali¬ 
brated using VIPERS data, while for stellar mass-limited sub- 
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Table 1. Definition of the samples. 


Redshift range 

Limiting magnitude 
M fl (z = 1.1) - 5 log(/i) 

N g 

W1 + W4 

(Tig 

S 3 g 

R = 8/T 1 Mpc 

S 4 g 

R = 8/z -1 Mpc 

0.5 <z 

< 0.1 

<-19.5 

8670 + 6863 

0.95 ± 0.06 

1.81+0.20 

8.13 + 2.03 

0.5 <z 

< 0.1 

<-20.0 

6101 +4963 

1.00 + 0.06 

1.82 + 0.22 

8.12 + 2.05 

0.5 <z 

< 0.1 

<-20.5 

3671 + 3025 

1.07 + 0.07 

1.83 + 0.27 

7.93 + 2.54 

0.5 < z 

< 0.1 

<-21.0 

1787 + 1478 

1.16 + 0.15 

1.78 + 0.30 

6.29 ± 2.96 

0.7 < z 

< 0.9 

<-20.0 

7455 + 5384 

1.01+0.05 

1.74 + 0.14 

7.19+1.34 

0.1 <z 

< 0.9 

<-20.5 

4979 + 3475 

1.05+0.04 

1.66 + 0.16 

6.02+1.38 

0.1 <z 

< 0.9 

<-21.0 

2457 + 1664 

1.10 + 0.06 

1.59 + 0.22 

5.50+1.76 

0.9 < z 

< 1.1 

<-20.5 

2751+1805 

1.12 + 0.07 

2.50 + 0.28 

14.11 ±3.10 

0.9 < z 

< 1.1 

<-21.0 

1752 +1067 

1.16 + 0.08 

2.54 + 0.38 

12.70 ±3.81 

Redshift range 

Limiting stellar mass 
log(M/M Q /T 2 ) 

N g 

W1 + W4 

(Tig 

S 3 g 

R = 8/z 1 Mpc 

S 4 g 

R = 8/z -1 Mpc 

0.5 <z 

< 0.7 

>9.0 

8745 + 6544 

0.97+0.10 

1.88 + 0.15 

8.51 ± 1.45 

0.5 <z 

< 0.7 

>9.5 

6091 +4318 

1.03+0.10 

1.94 + 0.15 

8.66+1.43 

0.5 <z 

< 0.7 

> 10.0 

3654 + 2581 

1.16 + 0.11 

2.02 + 0.16 

8.61 ± 1.45 

0.5 < z 

< 0.7 

> 10.5 

1292 + 713 

1.34 + 0.11 

1.90 + 0.18 

6.62+1.39 

0.7 < z 

< 0.9 

>9.5 

6159 + 4009 

1.09 + 0.08 

1.88 + 0.14 

7.59+1.37 

0.1 <z 

< 0.9 

> 10.0 

3746 + 2428 

1.18+0.08 

1.87 + 0.14 

7.29+1.36 

0.1 <z 

< 0.9 

> 10.5 

1467 + 819 

1.41 ±0.09 

2.04 ± 0.20 

7.81 ± 1.80 

0.9 < z 

< 1.1 

> 10.0 

1644 + 964 

1.23+0.08 

2.70 + 0.21 

13.28 ±2.43 

0.9 < z 

< 1.1 

> 10.5 

738 + 456 

1.43+0.09 

3.19 + 0.29 

16.18 ±3.88 


samples masses were assigned to galaxies us ing the stellar-to- 
halo mass relation (SHMR) of [Moster et al.l (l2013h . From the 
parent mock catalogues, a set of spectroscopic catalogues was 
derived by applying the same angular, photometric, and spectro¬ 
scopic selection functions as were applied to the real data. For a 
more detailed and complete description of the mock catalogues 
see lde la Torre et aD (120131) . 

>From the mock spectroscopic catalogues we derived 
volume-limited subsamples with cuts in blue absolute magni¬ 
tude and stellar mass corresponding to the observed ones. First 
of all, these mocks were used to test the effect of the gaps in the 
survey. As VIMOS is made of four quadrants 7x8 separated 
by 2 arcmin, characteristic cross-shaped gaps are left in the sur¬ 
vey; a further gap is present between the rows of pointings at 
different declination; finally, there are a few missing quadrants 
due to failed pointings. Cells whose projection on the sky in¬ 
cludes a gap can potentially miss some galaxies, which affects 
final counts. 

These gaps might be avoided by conservatively only count¬ 
ing galaxies in the cells that are completely included in one sin¬ 
gle quadrant, but in this way, only small scales would be sampled 
(the exact value obviously depends on the cell distance but it is 
generally lower than R ~ 5h~ l Mpc). Alternatively, the counts 
in each cell might be associated with the effective volume of the 
cell, subtracting the volume falling into the gaps; but this less 
drastic choice, which would slightly alter the shape of the cells, 
would still limit the range of the sampled scales. 

Another option would be filling the gaps. ICucciati et al.l 
(120141) applied two algorithms that use the photometric redshift 
information and assign redshifts to galaxies based upon the spec¬ 
troscopic redshifts of the nearest neighbours. In this way, it is 
also possible to take into account the varying completeness from 
field to field. Tests on mocks have shown that these algorithms 
are successful in reconstructing the lowest and highest density 
environments at a scale of 5 h~ l Mpc, but not in recovering the 
count PDF and its moments due to systematic biases. 



2 4 6 8 10 13 14 16 

R [h-> Mpc] 


Fig. 2. Fractional difference of the average £ 3 , £ 3 , S 3 and S 4 (from 
top to bottom) for the same set of mock catalogues as defined in Fig. 
1 , i.e. with 100 % sampling rate and without gaps, and with sampling 
rate and gaps as in VIPERS. The subsamples are limited at M B (z = 
1.1) - 51og(/r) < -20.5. Red triangles: 0.5 < z < 0.7; blue squares: 
0.7 < z < 0.9; green hexagons: 0.9 < z < 1.1. 


We therefore here adopted another solution. The tests on 
mocks have shown that when cells are not allowed to cross the 
gaps by more than 40% of their volume, the non-observed re¬ 
gions and the varying sampling rate can be approximated by a 
random Poisson sampling, and the original count PDF can be 
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Fig. 1. Comparison between mock catalogues with a sampling rate of 100% and without gaps (red triangles), and with a sampling rate and gaps as 
in VIPERS (blue triangles). The subsamples are limited at M B (z = 1.1) - 5 log(h) < -20.5. >From top to bottom: volume-averaged two-, three-, 
and four-point correlation functions, normalized skewness S 3 and kurtosis S 4 in redshift space. First column: 0.5 < z < 0.7; second column: 
0.7 < z < 0.9; third column: 0.9 < z < 1.1. 


recovered with good precision (Bel et al. 2015, in preparation). 
This means that to obtain good estimates of the quantities we dis¬ 
cuss here (/-point correlations and normalized moments), which 
depend on the density contrast Ap/p, it is sufficient to implement 
the restriction on the volume of the cells falling into the gaps. 

In our analysis, we conservatively only considered spherical 
cells for which no more than 30% of the volume falls in a gap. 
Moreover, to improve the statistics, we combined the counts of 
the W 1 and VV'4 fields. 

In Fig. |T] we show the results obtained from the analysis of 
mock subsamples limited at Mb(z - 1.1) - 5 log(/t) < -20.5 in 
the three redshift bins [0.5,0.7], [0.7,0.9], [0.9,1.1]. We compare 
the ideal case with 100% completeness and no gaps to the more 
realistic case with gaps and the same spectroscopic incomplete¬ 


ness as in our observed catalogue, that is, including the effects 
of the target sampling rate, TSR(Q), and the spectroscopic sam¬ 
pling rate, S S R(Q), where Q indicates the quadrant dependence. 

Two other selection effects were not taken into account: the 
colour sampling rate, CSR(z), and the small-scale bias due to 
the constraints in the spectroscopic target selection (slits cannot 
overlap). The first effect depends on redshift but it is weak in 
our redshift range (see Fig. 5 of lGuzzo et al.ll2014l) . while the 
second effect is negligible because the angular radii of our cells 
are generally larger than the size of one quadrant. 

We note that other sourc es of systematic errors, as discussed 
bv lHui & Gaztanagal (119991) , are the integral constraint bias, af¬ 
fecting the /-point correlation functions, and the ratio bias, af¬ 
fecting the estimate of S j. Given the large size of our volumes, 
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such systematic effects are weaker than the other errors, how¬ 
ever, and can be neglected. 

Figure Q] shows that the original values are recovered with 
good precision (within 1 cr error), particularly in the scale range 
between 4 and 10 h~ l Mpc. 

A more detailed analysis of the differences is possible with 
Fig.H which gives the fractional difference for ft, ft, ft, S 3, and 
4 as a function of scale for the same mock subsamples as in Fig. 
IH it shows that in most cases we can retrieve the 7-point correla¬ 
tion functions and S j with only a small systematic difference. In 
the first redshift bin (0.5 < z < 0.7) at a radius R = 8 h~ l Mpc, ft 
is overestimated by 8%, while ft is underestimated by 3% and 
ft by 6 %: this translates into an underestimate of 53 by 16% 
and of S 4 by 26%. We have similar values in the second red- 
shift bin (0.7 < z < 0.9). In the last redshift bin (0.9 < z < 1.1) 
the 7-point correlation functions show the largest difference, in¬ 
creasing with order 7: but these deviations at different orders are 
correlated, so that finally the values of S 3 at 8 h~ l Mpc is under¬ 
estimated by only 10 % and of 54 by 20 %, which is comparable 
to what is found for the other two redshift bins. The cause of 
the larger deviations in the last redshift bin is the lower density 
of the subsample; we take these systematics into account in the 
discussion of our results. 

It is interesting to point out that we find values between 
1.8 and 2.1 for .S '3 and between 8 and 10 for S 4 for mocks; 
as an example, the analysis of the mock subsamples limited at 
Mb(z = l.l)-51og(/z) < -20.5 in the redshift bin [0.7,0.9] gives 
S 3 ~ 2.13 + 0.16 and S 4 ~ 9.8 ± 1.6 at R = 8 /T 1 Mpc. S 3 and 
S 4 show no significant redshift evolution, and their values are 
also comparable within the errors to the value measured in local 
redshift surveys for galaxies in a similar luminosity range. 

Because we know both the cosmological and the “observed” 
redshift for galaxies in the mock samples, including the peculiar 
velocity and measurement error, we can estimate the conversion 
factor from redshift to real space from the mock samples. We 
need this factor to compare our results with second-order per¬ 
turbation theory predictions. Figure [3] shows the difference be¬ 
tween the estimates in real and redshift space for the subsamples 
limited at Mb(z = 1.1) - 5 log(/z) < -20.5 in the three redshift 
ranges. The redshift space correlation functions show the ex¬ 
pected loss of power at small scales and the reverse trend at large 
scales. The estimate of the volume-averaged two-point corre¬ 
lation function in redshift space is flatter than the correspond¬ 
ing estimate in real space; the difference becomes significant on 
scales smaller than ~ 4h 1 Mpc. While the real space values of 
S 3 and S 4 increase at smaller scales, the increase is suppressed 
in redshift space; the difference becomes small beyond ~ 4 h~ l 
Mpc. However, at small scales we have large errors due to the 
small number of objects in the cells. For these reasons we focus 
our analysis on the 4-10 h 1 Mpc range, and particularly at 8 
h 1 Mpc, where we expect to be in the quasi-linear regime and 
predictions of second-order perturbation theory should hold. 

We recall here another bias affecting mass-selected galaxy 
samples, which has_been discussed and tested with mock cata¬ 
logues bv lMarulli et al.l(l2013l) . The lowest stellar mass subsam¬ 
ples suffer from incompleteness because VIPERS is magnitude 
limited ( i AB < 22.5); as a consequence, we can miss high mass- 
to-ligh^atMgalaxi^^RmjLillS-SllSlSSiS-SLiiiS^l!-^^ 11 ^! 

et al. (120131 ? found that these galaxies are faint and red and that 
the clustering amplitude can be suppressed u p to 5 0% o n scales 
below 1 h~ l Mpc. However, as discussed bv lMarulli et al.l(l2013h . 
the abundance of red and faint galaxies is overpredicted by the 
semi-analytic model used for the tests, and the clustering of red 


galaxies appears to be overestimated with respect to real data (de 
la Torre et al. l201 lL Cucciati et al.l2012l) . so that the amplitude of 
the effect might be overestimated. As we have previously noted, 
we did not analyse small scales and did not correct for stellar 
mass incompleteness. 


5. Results 

5.1. Volume-averaged correlation functions 

In this section we present the results of our statistical analysis on 
the combined W 1 and W4 samples. 

Figure [4] shows the volume-averaged two-point correlation 
function obtained from counts in cells for luminosity- and stellar 
mass-limited subsamples in the three different redshift bins. 

In the same figure, as a reference for comparing the results 
in the different redshift bins, we plot the expected real space 
power-law ft in the redshift bin [0.5,0.7] for the M s (z = 1.1) - 
5 log(/z) < -20.5 subsamples (top panels) and M , > 1 0.0Aft 
(botto m panels), derived from the ft estimate of lMarulli et al.l 
d20l3h : we converted their two-point correlation function to the 
voIume-aver aged correlation function through the formula ( Pee- 
bles & Groth ll976lt : ~ 


ft =-ft. (17) 

52 2ft3 - y)(4 - y )(6 - yy~ 

The line shows the effects of redshift space distortions, which 
lower the value of ft on small scales and increase it on large 
scales. 

It is clear that the amplitude of ft increases with both lu¬ 
minosity and stellar mass at all redshifts. ft appears to have 
a stronger dependence on stell ar mass than on luminosity, in 
agreement with the results of lMarulli et al.l (! 2013 ): see their Fig.3 
for the redshift space two-point correlation functions. 

There are some fluctuations: for example, the dependence 
on luminosity appears to be sligthly weaker in the intermediate 
and distant redshift bins. However, these variations are consis¬ 
tent when taking into account statistical errors and sample vari¬ 
ance, which are included in error bars. We conclude that the de¬ 
pendence of the two-point correlation function on luminosity 
and stellar mass does not evolve significantly up to z ~ 1 . 

In Figs. [3 and [ 6 ] we show the volume-averaged three- and 
four-point correlation functions. Their behaviour reflects the 
two-point correlation functions, showing a stronger dependence 
of the correlation amplitude on stellar mass than on luminosity. 

The specific signature of the hierarchical scaling is the 
power-law relation between high-order correlation functions 
(Eq. U0l) . In Figs. [7] and [ 8 ] we show the three- and four-point 
volume-averaged correlation functions as a function of the two- 
point volume-averaged correlation functions. The data clearly 

follow the hierarchical scaling relations ft <x ft and ft oc ft. 
These relations appear to hold at all luminosities and masses in 
the the first two redshift bins, but some systematic differences 
appear in the last redshift bin, particularly for the stellar-mass 
limited subsamples, where points are systematically higher than 
the reference scaling law, but in this case the values are also con¬ 
sistent with the same scaling relation observed at lower redshifts. 

As we have previously discussed, the existence of these scal¬ 
ing relations has been verified in the local Universe: they are ex¬ 
pected for the matter distribution in the quasi-linear regime, as 
a consequence of gravitational clustering. In this case, it is nat¬ 
ural that they do not evolve with redshift: however, it is not an 
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Fig. 3. Comparison between mock catalogues with 100% sampling rate and without gaps in real space (red squares) and redshift space (blue 
triangles). The subsamples are limited at M B {z = 1.1) —5 log(/t) < -20.5. First column: 0.5 < z < 0.7; second column: 0.7 < z < 0.9; third column: 
0.9 < z < 1.1. 


obvious result to observe the same hierarchical behaviour for the 
galaxy distribution at all redshifts, given the evolution of bias. 


5.2. Skewness and kurtosis 

>From the counts in cells we derived the rms cr (Eq.[3}, the nor¬ 
malized skewness S 3 and kurtosis S 4 (Eq. ITU for the different 
VIPERS subsamples. Their values at R = 8 h~ l Mpc are given in 
Cols. (4), (5), and ( 6 ) of Table [T] The R = 8 /C 1 Mpc reference 
radius is nearly optimal because it is large enough to enter into 


the quasi-linear regime, and at the same time it is in the scale 
range for which we have a good sampling. 

In Figs. |9]and[T0]we show S 3 and .S ’4 as a function of lumi¬ 
nosity and stellar mass in the three redshift bins. We also show 
the predictions of second-order perturbation theory in real space 
for the matter distribution and the corresponding predictions for 
galaxies, derived from the matter value assuming the linear bias 
estimated from £ 2 , and corrected for redshift space distortion us¬ 
ing the factors obtained from mocks. This derivation is described 
in the next subsection. The theoretical curves for S 3 and .S’ 4 are 
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s [h _1 Mpc] s [hr 1 Mpc] s [h -1 Mpc] 



Fig. 4. Volume-averaged two-point correlation functions £ 2 as a function of the B absolute luminosity (upper panels) and stellar mass (lower 
panels). The limits of the subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels. Black triangles: 
M B {z = 1.1) — 5 log(/7) < -19.5 (log(M»/M 0 /t -2 > 9.0A/ G ); red squares: M B (z = 1.1) — 5 log(/j) < -20.0 (log (M,/M 0 hr 2 > 9.5M 0 ); blue pentagons: 
M b {z = 1.1) — 5 log(/r) < -20.5 ( log(M,/M e lr 2 > 10.0A/ o ); green hexagons: M B (z = 1.1) - 5 log (h) < -21.0 (log {M,/M Q h~ 2 > 10.5A/ o ); magenta 
heptagons: M B (z = 1.1) - 51og(/i) < -21.5 (log (M,/M e h~ 2 > ll.OM 0 ). Dotted lines: real-space for the M B (z = 1. 1) ^51og(/i) < -20.5 
(log(M,/M 0 /r 2 > 10.0M o ) subsamples in the redshift bin [0.5,0.7], predicted from the power-law fit of f 2 in lMarulli et ali (2013 ). 


shown for radii larger than ~ 6 h 1 Mpc, as they are calculated in 
the quasi-linear regime. 

In the first redshift bin, both for luminosity and stellar mass 
limited samples, the value of S 3 is constant and around 2 at 
small and intermediate scales, but it starts decreasing beyond 
R ~ 8 h~ l Mpc. In principle, variations of .S '3 with scale can 
be due to changes in the slope of the power spectrum or to a 
scale-dependent bias. However, such a systematic effect can be 
ascribed to the small number of independent cells at large scales, 
as shown by mocks and reflected in the large error bars. In the 
same redshift bin, .S '4 shows a small decrease at large scales and 
is consistent with a constant value of ~ 7.3 between 4 and 10 


hr 1 Mpc. In the range 6 - 8 !r x Mpc, the best scales to com¬ 
pare with perturbation theory (on larger scales the errors in¬ 
crease significantly), the theoretical predictions for S 3 and .S '4 
are slightly higher than the observed values corresponding to the 
Mb(z = 1.1) - 5 log(/r) < -20.5 subsample, but only at l<x level. 

In the second redshift bin the value of S 3 for luminosity- 
limited subsamples is around 1 . 8 , sligthly lower than in the first 
bin, but still consistent within the errors; moreover, it is consis¬ 
tent with a constant value in the whole range of scales. The value 
of S 3 for mass-limited subsamples is also constant in the whole 
range of scales and is consistent with the value in the first red¬ 
shift bin. S 4 has an analogous behaviour: while showing a sys- 
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Fig. 5. Volume-averaged three-point correlation functions as a function of the B absolute luminosity (upper panels) and stellar mass (lower panels). 
The limits of the subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels (symbols and colours are the 
same as in Fig.0. 


tematic decrease, particularly in luminosity-limited subsamples, 
it is still consistent with a constant value in the range 4 - 16/z _1 
Mpc. As in the case of the first redshift bin, in the range 6 - 8 /z -1 
Mpc the theoretical predictions for S 3 and 54 are slightly higher 
than the corresponding observed values. 

In the third redshift bin the values of S 3 and S 4 for 
luminosity- and stellar-mass-limited subsamples increase sys¬ 
tematically with scale. Moreover, in contrast with the two pre¬ 
vious redshift bins, in the range 6 - 8 /z~' Mpc, the theoretical 
predictions for S 3 and .S '4 are lower than the observed values. 

To better appreciate the significance of these deviations, we 
note that of 26 mocks, 3 show an increase of the values of S 3 
and S 4 similar to what we find in the last redshift bin. 

In fact, higher-order statistics are very sensitive to large- 
scale structure, and the correlated variations in the measured 


values of S 3 and .S '4 probably indicate genuine fl uctuati ons in 
the galaxy distribution (see e.g. the discussion in lCroton et al.l 
l2004bl) . 

In our case, this interpretation is suggested by checking the 
W 1 and W4 fields separately: we find that in the outermost red¬ 
shift shell, both S 3 and S 4 are larger in W 1 than in VV’4. For 
example, for the Mb(z = 1.1) - 5 log(/z) < -21.0 subsample, at 
R =8 /T 1 Mpc, we find S 3 = 2.7 ± 0.5 in IV1 and S 3 = 1.6 ± 0.3 
in W4-. Analogously, for the log {M/M Q h~ 2 ) > 10.5 subsample 
at R =8 li -1 Mpc, we find S 3 = 3.4 + 0.5 in the W 1 field and 
S 3 = 2.0 ±0.3 in the W4 field. This difference might be regarded 
as the imprint of spatially coherent structures more prominent in 
Wl. 

In conclusion, the values of S 3 and S 4 do not show any sig¬ 
nificant dependence on luminosity or on stellar mass: the points 
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Fig. 6. Volume-averaged four-point correlation functions as a function of the B absolute luminosity (upper panels) and stellar mass (lower panels). 
The limits of the subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels (symbols and colours are the 
same as in Fig. 0- 


corresponding to different subsamples are consistent within the 
error bars (we discuss a possible weak dependence on luminos¬ 
ity in the next subsection). There is no evidence of evolution in 
redshift either, apart from the systematic increase of S 3 and S 4 
with scale in the last redshift bin. 

Taking into account the behaviour of mocks, the observed 
systematic variations in the values of high-order moments are 
consistent with the fluctuations expected for comparable vol¬ 
umes randomly extracted from a ACDM universe. 

It is possible to compare our results on 5 3 and 5 4 with those 
obtained bv lWolk et al.l (12013l) for the four CFHTLS-Wide fields. 
They have divided the galaxies in the photometric catalogue into 
four redshift bins through the estimated photometric redshifts; 
for galaxies with M g < -20.7, they have estimated S j as a func¬ 
tion of angular scale and the corresponding 3D values through 


deprojection, which, as they discussed, rely on some approxima¬ 
tions. Their work is therefore complementary to ours: they have 
a larger area and number of objects, but we can directly estimate 
the 3D (redshift space) S j; they can sample smaller, highly non¬ 
linear scales where we do not have enough statistics, but we can 
better sample the quasi-linear scales; finally, we can also test the 
dependence of S 3 on luminosity and stellar mass. 

A comparison with their Fig. 12 shows that, as expected (see 
our Fig[3}. their deprojected values for S 3 and .S '4 on small scales 
(R <5 h~' Mpc) are higher than our redshift space values. On 
larger scales, the redshift space effect on S 3 and .S '4 becomes 
negligible, and their estimate is consistent with ours. 

We note that lWolk et al.l (12013l) found significant deviations 
in the results for the IV3 field, while we have found differences 
between 1V1 and 1V4 in our last redshift bin: this shows that sam- 
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Fig. 7. Scaling relation of the volume-averaged two- and three-point correlations as a function of the B absolute luminosity (upper panels) and 
stellar mass (lower panels). The limits of the subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels 
(symbols and colours are the same as in Fig. [4}. The dashed line represents the scaling relation log 10 (£ 3 ) = 2 log 10 (£ 2 ) + log 10 (2). 


pie variance is still significant for high-order statistics on the 
scale of CFHTLS Wide Fields. 


5.3. Implications for biasing 

We now discuss the implications of our analysis for biasing. We 
concentrate on the reference scale R = 8/z 1 Mpc, where second- 
order perturbation theory predictions can be applied and results 
are still reliable (errors and systematic deviations increase on 
larger scales). Because we aim to compare our results with the 
matter density field, statistical quantities referring to galaxies are 
indicated with a subscript g and those relative to matter with a 
subscript m. 


Figure fill shows the values of cr$ g (top panel) and S z g (bot¬ 
tom panel) at R — 8/z 1 Mpc for the VIPERS volume-limited 
subsamples with different limiting absolute magnitudes and in 
the different redshift bins. In the same figure we also show the 
corresponding VVDS estimates (IMarinoni et al . 2005) and the 
2dFGRS estimates for the local Universe dCroton et al.l 12004^ ) 
for galaxies with a similar luminosity as ours. 

At a given redshift, VIPERS subsamples with a brighter ab¬ 
solute magnitude limit have higher values of cr% g , but there is no 
significant evolution of cr& g with redshift. The same holds when 
combining our results with those of the 2dFGRS in the local 
Universe and those of the VVDS at higher redshift: cr& g shows 
no significant evolution from z - 0 to z = 1.4 (VVDS points 
are systematically lower but at the lcr level). This implies (see 
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Fig. 8. Scaling relation of the volume-averaged two and four-point correlations as a function of the B absolute luminosity (upper panels) and 
stellar mass (lower panels). The limits of the subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels 
(symbols and colours are the same as in Fig. 0. The dashed line represents the scaling relation log 10 (£ 3 ) = 3 log lofe) + l°gio( 8 )- 


e.g. the discussion in lMarinoni et al.f2005l) a strong evolution of 
the linear bias b with redshift because cr 8m increases with time 
(see Eq. O There are various models that describe the evolu¬ 
tion of^^)and i ex 2 l 3 Mijfs i< decrease > withJime_(see_e i ^.Blanton 
et al. l 2000 h : from an empirical point of view, we note that the 
available data can be fitted by the simple relation b(z) or 1 /cr 8m . 

The skewness S i g of the VIPERS subsamples measured at 
8 h~ l Mpc and plotted as a function of redshift has more fluc¬ 
tuations than crsg, with a minimum value in the redshift bin 
[0.7,0.9], but it does not show a significant dependence on lumi¬ 
nosity and is still consistent with a constant value independent of 
redshift. The values of S i g in the V VDS below z = 1.2 are lower 
than VIPERS values, but are consistent within the errors, while 
they start to decrease beyond z ~ 1 . 1 . 


The absence of a significant evolution of S 3 g with redshift 
is not limited to our redshift range: the values of S 3 g measured 
in VIPERS are similar to those measured in the 2dFGRS, that 
is, S 3 ~ 2.0 ± 0 . 2 , where depending on the subsample S 
varies from 1.95 to 2.58 (while not shown in the figure, the val¬ 
ues of S 4 are also consistent with the 2dFGRS ones). There¬ 
fore, taking into account all data points, starting from the lo¬ 
cal value for the 2dFGRS up to z = 1.1 (VIPERS and VVDS 
data), S 3 g is consistent with a constant value ~ 2: in VIPERS 
the strongest but marginal deviations of the S 3 g value are for 
Mg(z - 1.1) - 5 log(/t) < -20.0 galaxies in the nearest redshift 
range [0.5,0.7] and for Mg(z = 1.1) — 5 log(/j) < -21.0 galaxies 
in the most distant redshift interval [0.9,1.1], both giving a value 
of S 3 g that is 15% higher. 
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s [h _1 Mpc] s [hr 1 Mpc] s [h -1 Mpc] 



Fig. 9. Normalized skewness S 3 as a function of the B absolute luminosity (upper panels) and stellar mass (lower panels). The limits of the 
subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels (symbols and colours are the same as in Fig. [4] We 
note that black points corresponding to galaxies with M < -19.50 or M, > 9.0Af 0 are only plotted for the first redshift bin ([0.5, 0.7]), but most 
of them are not visible as they lie below the points of the other samples. The dashed line corresponds to S 3 = 2. In the top panels, the dotted line 
is the prediction from second-order perturbation theory for the matter distribution in real space; the solid line is the prediction of 53 for galaxies 
with M b (z = 1.1) — 51og(/r) < -20.5 (to be compared to blue pentagons): it was obtained from the matter values, converted to redshift space and 
divided by the corresponding linear bias factor. 


Figurefl2lshows crg 9 (top panel) and S 3 g at 8 /z 'Mpc (bottom 
panel) as a function of absolute magnitude for the three redshift 
bins. cr% g shows a systematic increase with luminosity (reflecting 
the dependence of the correlation amplitude on luminosity), but 
at a given absolute luminosity its value is similar in the three 
redshift bins. 

S 3 g appears to be independent of absolute magnitude, with 
fluctuations from sample to sample. However, if we exclude the 
points relative to the last redshift bin, where S 3 g has a higher 
value, the data might suggest a small decrease of S :\ g with in¬ 


creasing luminosity, reminiscent of the results of lCroton et al.l 
(l2004bf) for the 2dFGRS. 

A trend of S 3 g with luminosity is interesting because in the 
hypothesis of linear biasing, S i g is inversely proportional to the 
bias factor b: knowing from the two-point correlation function 
of our samples that b increases with luminosity, we expect a cor¬ 
responding decrease of Si g . 

To test whether our results are consistent with the linearity 
of bias, we therefore estimated the bias of galaxies with respect 
to the underlying matter density field at R - 8 /z 1 Mpc, using 
the observed crs g and S j, g of the galaxy distribution and estimat- 
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s [h _1 Mpc] s [hr 1 Mpc] s [h -1 Mpc] 



Fig. 10 . Normalized skewness S 4 as a function of the B absolute luminosity (upper panels) and stellar mass (lower panels). The limits of the 
subsamples in absolute magnitude and stellar mass are shown in the left upper and lower panels (symbols and colours are the same as in Fig. 0 . 
We note that black points corresponding to galaxies with M < -19.50 or M, > 9.0 M e are only plotted for the first redshift bin ([0.5, 0.7]), but most 
of them are not visible as they lie below the points of the other samples. The dashed line corresponds to £4 = 8 . In the top panels, the dotted line 
is the prediction from second-order perturbation theory for the matter distribution in real space; the solid line is the prediction of S 4 for galaxies 
with M b (z = 1.1) - 5 log(/i) < -20.5 (to be compared to blue pentagons): it was obtained from the matter values, converted to redshift space and 
divided by the corresponding linear bias factor. 


ing cr m and S 3m of the matter distribution through perturbation 
theory. 

Juszkiewicz et al. (1 19931) and lBernardeaul (1 1 994allbh (see also 
iBernardeau et al.ll2002l and references therein) have shown that 
for a smoothed density field with primordial Ga ussi an fluctua¬ 
tions, Peebles’ unsmoothed value of S 3m = 34/7 dPeebleJl980h 
has to be corrected according to the expression 


S 3m = 34/7 + d In crl/d In R, (18) 


where r/ln o-^/dlnR is the logarithmic slope of the linear 
variance of the matter density field smoothed with a spherical 
top-hat function of radius R, 

1 r° 

<T 2 m {R) = — 2 dkk 2 P(k)W 2 (kR). (19) 

For a power-law spectrum P(k ) oc k", Ea.ll 8 lbecomes S 3m = 
34/7 - (n + 3). 

Similar relations hold for higher orders, involving higher- 
order derivatives. 
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Fig. 11. o ~(upper panel) and S 3g (lower panel) as a function of red- 
shift. Colours represent the different absolute magnitude limits as in pre¬ 
vious figures. The open hexagon represents the corresponding value for 
2dFGRS ealaxiesjn = &e ii absolutejna2nitude i Mj u _ran2e_[ ; 21 j; 20 i lJ|Cro- 
ton et al. l2004bfr Black op en squares are the va lues for VVDS galaxies 
brighter than M B = -21 dMarinoni et al.l[20051) . For better visibility, 
points corresponding to different redshifts are slightly shifted in magni¬ 
tude. 


The values obtained from perturbation theory have been 
tested with numerical simulations, and it has been shown that 
in the range we are studying, that is, at R = 8/T 1 Mpc and for 
cr 8m ~ 1 , they are very accurate: for example, the difference in 
the S 3 values is smaller than a few percent ([ Baugh e t al.ll 19951 

iFosalba & Gaztanagalll998LlBernardeau et al.ll2002[). 

Applying Eqs. (fTST l and (119b and using the software CAMB 
jLewis & Bridlel 12002 ). we have computed the values of cr 8m 
and S 3m for a power spectrum with the new cosmological pa- 
ram eters deri ved fro m the Planck mission ( Planck Collaboration 
et al. teOOal) ' and with the old Millennium parameters (first year 
WMAP data and 2dFGRS, with O m = 0,25, Q A = 075, n = 1 
and o- 8m = 0.9). 


We here assumed that the standard KCDM model is correct. 
With other assumptions, such as a dark energy component with 
an evolving equation of state or modified gravity, the cl uste r¬ 
ing an d bias evolution would be affected (see e.g. M unshi e t alJ 
:2004 ). as would the redshift distortions (iHellwing et al.ll2013h . 
This dependence on cosmology will be studied in a future work. 

We also converted the observed cr% g and S 3g to real space 
values by applying correction factors directly derived from the 
mocks. 

For the subsample limited at Mg(z — 1.1) — 5 log(/z) < 
-20.50, we give in Table 2 the redshift range (Col. 1), the values 
of cr & g (Col. 2), cr 8m (Col. 3), b = (r ig /cr m (Col. 4), S 3g (Col. 
5), S 3 ,„ (Col. 6 ), all measured at a scale of R - 8/T 1 Mpc. 

In Fig. [13] we plot our estimates for the linear bias term b 
as a function of redshift, with the c orreponding es timates for 
VIPERS of iMarulli et ahl (floTl) and lDi Porto etaD (|20I4|) . As 
expected, the estimates are fully consistent, with b in creas ing 
with luminosity and redshift. As discussed by IDi Porto et all 


Fig. 12. cr Sg (upper panel) and S 3g at R=8 h 1 Mpc (lower panel) as a 
function of galaxy luminosity for the three redshift bins. Red triangles: 
0.5 < z < 0.7; blue squares: 0.7 < z < 0.9; green hexagons: 0.9 < z < 
1.1. For better visibility, points corresponding to different redshifts are 
slightly shifted in magnitude. 


(120141) . there is only a difference in the last redshift bin where 
the estimate of Marulli et al. is lower than that of Di Porto et al.l 
(1201 4) . The difference is probab ly due to the way b i s estimated 
(c ounts in cells in our c ase and in lDi Porto et al.l2014l pair counts 
in IMarulli et al Jl 2 () 1 3l) . Our estimate is consistent with both the 
other two estimates at the lcr level, however. 

In Fig.[l4]we compare the linear bias directly measured from 
the ratio of the galaxy and matter rms, b = cr 3g lcr 3m , with the 
ratio of the galaxy and matter skewness, S 3 m /S 3g . Under the hy¬ 
pothesis of linear biasing, the two ratios should have the same 
value. For the first two redshift bins we find slightly different val¬ 
ues: the skewness ratio is systematically higher than the bias di¬ 
rectly computed from the variance. The third redshift bin shows 
the largest discrepancy, but with the opposite behaviour, that is, 
the skewness ratio is lower than the bias directly computed from 
the variance. This different behaviour is a consequence of the 
fact that the value of S 3g in the last redshift bin increases with 
scale and becomes higher than at lower redshifts. 

We can quantify the degree of non-linearity by directly esti¬ 
mating the second-order term b 2 from Eg. 1161 


b 2 = KbS 3g -S 3m )/3, (20) 

where we used the real space values S 3g and b obtained from 
the redshift space values by using the conversion factor calcu¬ 
lated from the mocks. We note that this correction is small (a 
few percent) at our scale of R - 8 h~ l Mpc, because this scale is 
at the transition from the regime of small-scale velocity disper¬ 
sion (where redshift space correlation functions are lower than 
real space ones) to the regime of infall where redshift space cor¬ 
relations are higher than real space ones (see Fig. [3} 

In this formalism, if b > 0, b 2 is negative when crg g /cr 3m < 
S 3m /S 3g . This is what happens in the first two redshift bins, 
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Table 2. Real space values of cr 8 and S 3 of galaxies with M B (z = 1.1) - 51og(/i) < -20.50 vs. those expected from second-order perturbation 
theory (with Millennium and Planck cosmological parameters). 


Redshift range 

Cig 

(real space) 

O' 8/7? 

(WMAP/Planck) 

b 

(linear bias) 

Si g (R - 8/z 1 Mpc) 
(real space) 

S 3 m (R = 8 h 1 Mpc) 
(WMAP/Planck) 

0.5 < z < 0.7 

1.00 + 0.12 

0.61/0.63 

1.65 + 0.20 

1.91 ±0.29 

3.52/3.48 

0.7 < z < 0.9 

0.98 ± 0.08 

0.55/0.57 

1.77 + 0.14 

1.76 + 0.18 

3.52/3.48 

0.9 <z < 1.1 

1.04 + 0.13 

0.51/0.52 

2.04 ± 0.27 

2.28 ± 0.25 

3.52/3.48 



0.4 0.6 0.8 1 1.2 

Redshift 



M b (z= 1.1) - 51og(h) 


Fig. 13. Linear bias b as a function of redshift. Top panel: M B < 
-20.0(z = l.l) + 5 log(/r); middle panel: M B < -20.5(z = 1.1) + 5 log(/r); 
bottom panel: M B < - 21 ,0(z = 1.1) + 51og(/i). Red he xagons: ouresti- 
mates of b = <x R „/< x R ,„. Blue squares: est imates of lMarulli et all d2013t) . 
Magenta triangles: iDi Porto et al .(Mi- 


Fig. 14. Estimates of b = <Tg g /crg m (solid lines, filled symbols) and b' = 
Sim/Sig (dotted lines, open symbols); b’ = b under the assumption of 
linear bias. Red lines with triangles: 0.5 < z < 0.7; blue lines with 
squares: 0.7 < z < 0.9; green lines with hexagons: 0.9 < z < 1.1. 


where at nearly all magnitudes £>2 is negative: for example, for 
the subsample limited at M < -20.5(z=l.l) - 5 log(h), we find 
(?2 = -0.20 ± 0.49 in the first redshift bin and (72 = -0.24 ± 0.35, 
in the second redshift bin. In contrast, we find a positive bi in the 
third bin, with £>2 = +0.78 ± 0.82. 

As we have noted above when discussing the results of our 
tests on mocks, the assumption that masked regions and inho¬ 
mogeneities can be described as a Poissonian random sampling 
gives a small bias with an overestimate of b of a few percent and 
an underestimate of S3 around 10-15%. Using the correction fac¬ 
tors derived from the average of the mocks, we find for the sub¬ 
sample limited at M < -20.5(z=l.l)-5 log(h)(?2 = -0.03+0.49 
in the first redshift, bi = -0.25 ± 0.35 in the second redshift bin, 
and f>2 = +0.72 + 0.82 in the third bin. The differences are well 
within 1 cr error. 

It would be tempting to interpret these results as suggest¬ 
ing a possible evolution of the non-linear bias Z?2 with redshift, 
with a simi lar trend, for example, as for the model of Sefusatti 
& Komatsu (120071) . Unfortunately, the problem is the extreme 
sensitivity of £>2 to the errors on b and S j, g , amplified by a factor 
b 2 , and we have seen that subsamples in the last redshift bin are 
affected by larger errors and systematic trends. 

With these caveats, we can check the consistency of our 
results with other works in the same redshift and luminosity 


ranges. In this comparison, one has to take into account the sen¬ 
sitivity o f bi to the different methods and, as pointed out by Ko- 
vac et al. 1201 1[), to sample variance. In fact, even local measure- 
ments of the non- lin ear term have given dif ferent values (see e.g. 
Verde et al.ll2002l andlPan & Szapudill2005l and the discussion in 


Gaztanaga et akll2005^ndjMarinoni et al.ll2008l). 


First of all, Di Porto et al.1 ( 20141) have analysed the VIPERS 
data reconstructing the bias relation from the estimate of the 
probability distribution function: they found a small (< 3%) but 

significant deviation from linear bias. _ 

In their analysis of the four CFHTLS Wide fields. fWolk et al.l 
( 2013 !) have found that perturbation theory predictions agree 
well with their measurements when taking into account the lin¬ 
ear bias, but note that there is still a small discrepancy that can 
be explained by the presence of a non-linear bias term. This is 
alsoconsistent with what we found. 

Marinoni et al. (2005 ) have analysed VVDS volume-limited 
samples limited at Mb < -20 + 5 log(/z) (this limit was fixed and 
did not take into account luminosity evolution) in the redshift 
bins 0.7 < z < 0.9 and 0.9 < z < 1.1, finding £>2 = -0.20 + 0.08 
and bi = -0.12 ± 0.08 (here the errors do not include sample 
varia nce): these v alues are consistent with ours below z = 1 . 

Kovac et al. (1201 ll) analysed the zCOSMOS galaxy overden¬ 
sity field and estimated the mean biasing function between the 
galaxy and matter density fields and its second moment, finding 
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a small non-linearity, with the nonlinearity parameter b/b (de¬ 
fined in the formalism of lDekel & Lahavll 19991) at most 2% with 
an uncertainty of the same order. 

Gaztanaga et al. (l2005h have found b = 0.93^'qj and C 2 = 
b 2 /b — -0.34^-Igfrom the measurement of the Q 2 parameter in 
the three-point correlation function of the 2dFGRS for the local 
Universe. 

The non-linear term we have measured in the redshift inter¬ 
val between z = 0.5 and z — 0.9 is therefore similar to what has 
been measured in the above surveys. We conclude that there is 
general evidence for a small but non-zero non-linear b 2 term. It 
is also clear that no evolution of hi with redshift can be detected 
in the available data, in contrast to the linear bias term. 

6. Conclusions 

We have analysed the high-order clustering of galaxies in the 
first release of VIPERS, using counts in cells to derive the 
volume-averaged correlation functions and normalized skew¬ 
ness S 2g and kurtosis S 4 g . We have analysed volume-limited 
subsamples with different cuts in absolute magnitude and stel- 
la r mass i n three redshift bins; these subsamples are the same as 
in lMarulli et all J20l3h . 

Errors were estimated through a set of mock catalogues, 
derived from dark mat ter halo cat alogu es repopulated with the 
method of Ide la Torre & Peacockl (i2013f) . The mocks were built 
to reproduce the properties of VIPERS, including masks and se¬ 
lection effects. Our analysis has shown that the high-order statis¬ 
tical properties of these mocks are consistent with observations. 

We also studied the dependence of the second- and third- 
order statistics of galaxy counts on the bias, deriving the linear 
bias term b and the first non-linear term hi, and comparing our 
results with predictions from perturbation theory and with other 
works in the literature. 

Here are our main conclusions. 

- We showed that the hierarchical scaling relations oc and 

- -3 

£4 cc £ 2 hold in the range of scales and redshifts we could 
sample, that is, 3 < R < 10/z 1 Mpc and 0.5 < z < 1.1. These 
relations are consistent with predictions from gravitational 
clustering and with the scaling observed in local surveys. 

- S 2g and S 4 g appear to be independent of luminosity; how¬ 
ever, if we do not not take the last redshift bin into account, 
there is a slight decrease of S 2g with increasing luminosity, 
an effectnreviousl^_detected > locall^jii i the < 2dFGRS i b^ i Cro¬ 
ton et al. (I2004bl) . 

- The values of S 2g and S 4 g are scale-independent within the 
errors and do not evolve significantly at least up to z = 0.9. 
We detected a systematic increase with scale in the last red¬ 
shift bin (beyond ~ 10/; -1 Mpc), mainly due to one of the 
two CFHTLS fields (Wl); this deviation is consistent with 
what can be expected from the sample variance shown by 
mock catalogues. 

- The observed values of ^3^ ~ 2 ± 0.2 and S 4 g ~ 8 + 0.4 
are similar to those measured in local surveys for galaxies 
in the same luminosity range. This confirms the substantial 
absence of evolution of S 2g in the redshift range 0 < z < 1 
at the level of ~ 10%. This result is expected for S 3 ,,,, but is 
not trivial for S 2g , given the evolution of bias. 

- At second order, galaxies with higher luminosity or stellar 
mass have a larger amplitude (greater linear bias parame¬ 
ter) of the volume-averaged two-point correlation function, 
consistently with the direct analysis of the two-point cor¬ 
relation function by [Marulli et al. ( 2013 ). We showed that 


our estimate of the linear bias parameter b - cr^J cr X) „ is 
con sistent within l c r with those of iMarulli et al.l (1201 3h 
and iDi Porto et al.l (120141) . The linear bias increases both 
with luminosity and with redshift: in our redshift range, 
we measured the lowest bias b — 1.47 + 0.18 for galaxies 
with Mb(z - 1.1) - 51og(/i) < -19.5 in the redshift bin 
0.5 < z < 0.7 and the largest bias b - 2.12 ± 0.28 for galax¬ 
ies with Mb(z = 1.1) — 5 log(/r) < -21.0 in the redshift bin 
0.9 <z< 1 .1. 

- For a given luminosity class, a% g does not evolve with red¬ 
shift. For example, comparing our values for Mb(z - 1.1) - 
5 log(/i) < -20.5 to the corresponding value measured in 
the 2dFGRS, we found that crg g is consistent with a constant 
value 1.0 (our lcr error is 10%), from z = 0 to z ~ 1. Given 
that cr 8m increases with time, we have the empirical relation 
b(z) oc l/cr 8m (z). 

— The value of the non-linear bias parameter b 2 measured be¬ 
low z ~ 1 at the scale R - 8 hr 1 Mpc, that is, in the quasi- 
linear regime, is negative but not statistically different from 
zero when taking into account the error; however, taking into 
account the ensemble of results coming from this_and other 


surveys in the redshift range 0.5 < z < 1 dMarinoni et al. 

2005 

. iKovac et al.ll201lL IWolk et alJ 2013. Di Porto et al. 

2014 

term 

evoli 

data. 

), there is evidence for a small but non-zero non-linear 
Including the results from local surveys as well, no 
ition of b 2 with redshift can be detected in the available 


- The comparison with the properties of mocks and with the 
predictions of perturbation theory shows that our results are 
consistent with the general scenario of biased galaxy for¬ 
mation and gravitational clustering evolution in a standard 
AC DM cosmology. 

In conclusion, we have provided an independent check on 
the second-order statistical studies of the galaxy distribution 
through our analysis; we explored the galaxy bias with an inde¬ 
pendent technique; finally, we determined the higher-order sta¬ 
tistical properties of the galaxy distribution in the redshift range 
between 0.5 and 1.1, thanks to the combination of volume and 
density of galaxies in the VIPERS survey. When VIPERS is 
complete, it will be possible to perform a more general analysis, 
which will allow us not only to decrease error bars, but also to 
include the dependence of high-order statistics on galaxy colour, 
to apply other high-order statistical tools such as the void prob¬ 
ability function, and to give better constraints on the non-linear 
bias. 
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